SNANA: A Public Software Package for Supernova Analysis

We describe a general analysis package for supernova (SN) light curves, called SNANA, that contains a simulation, light curve fitter, and cosmology fitter. The software is designed with the primary goal of using SNe Ia as distance indicators for the determination of cosmological parameters, but it can also be used to study efficiencies for analyses of SN rates, estimate contamination from non-Ia SNe, and optimize future surveys. Several SN models are available within the same software architecture, allowing technical features such as K-corrections to be consistently used among multiple models, and thus making it easier to make detailed comparisons between models. New and improved light-curve models can be easily added. The software works with arbitrary surveys and telescopes and has already been used by several collaborations, leading to more robust and easy-to-use code. This software is not intended as a final product release, but rather it is designed to undergo continual improvements from the community as more is learned about SNe. Below we give an overview of the SNANA capabilities, as well as some of its limitations. Interested users can find software downloads and more detailed information from the manuals at http://www.sdss.org/supernova/SNANA.html .Comment: Accepted for publication in PAS

The Core Collapse Supernova Rate from the SDSS-II Supernova Survey

We use the Sloan Digital Sky Survey II Supernova Survey (SDSS-II SNS) data to measure the volumetric core collapse supernova (CCSN) rate in the redshift range (0.03<z<0.09). Using a sample of 89 CCSN we find a volume-averaged rate of (1.06 +/- 0.19) x 10**(-4)/(yr Mpc**3) at a mean redshift of 0.072 +/- 0.009. We measure the CCSN luminosity function from the data and consider the implications on the star formation history.Comment: Minor corrections to references and affiliations to conform with published versio

Improving the LSST dithering pattern and cadence for dark energy studies

The Large Synoptic Survey Telescope (LSST) will explore the entire southern sky over 10 years starting in 2022 with unprecedented depth and time sampling in six filters, $ugrizy$. Artificial power on the scale of the 3.5 deg LSST field-of-view will contaminate measurements of baryonic acoustic oscillations (BAO), which fall at the same angular scale at redshift $z \sim 1$. Using the HEALPix framework, we demonstrate the impact of an "un-dithered" survey, in which $17\%$ of each LSST field-of-view is overlapped by neighboring observations, generating a honeycomb pattern of strongly varying survey depth and significant artificial power on BAO angular scales. We find that adopting large dithers (i.e., telescope pointing offsets) of amplitude close to the LSST field-of-view radius reduces artificial structure in the galaxy distribution by a factor of $\sim$10. We propose an observing strategy utilizing large dithers within the main survey and minimal dithers for the LSST Deep Drilling Fields. We show that applying various magnitude cutoffs can further increase survey uniformity. We find that a magnitude cut of $r < 27.3$ removes significant spurious power from the angular power spectrum with a minimal reduction in the total number of observed galaxies over the ten-year LSST run. We also determine the effectiveness of the observing strategy for Type Ia SNe and predict that the main survey will contribute $\sim$100,000 Type Ia SNe. We propose a concentrated survey where LSST observes one-third of its main survey area each year, increasing the number of main survey Type Ia SNe by a factor of $\sim$1.5, while still enabling the successful pursuit of other science drivers.Comment: 9 pages, 6 figures, published in SPIE proceedings; corrected typo in equation

Enhanced CP Violation with B→KD0(D‾0)B\to K D^0 (\overline D^0) Modes and Extraction of the CKM Angle gamma

The Gronau-London-Wyler (GLW) method extracts the CKM angle $\gamma$ by measuring $B^\pm$ decay rates involving $D^0/\overline D^0$ mesons. Since that method necessitates the interference between two amplitudes that are significantly different in magnitude, the resulting asymmetries tend to be small. CP violation can be greatly enhanced for decays to final states that are common to both D^0 and $\overline D^0$ and that are not CP eigenstates. In particular, large asymmetries are possible for final states f such that $D^0\to f$ is doubly Cabibbo suppressed while $\overline D^0\to f$ is Cabibbo allowed. The measurement of interference effects in two such modes allows the extraction of $\gamma$ without prior knowledge of $Br(B^-\to K^- \overline D^0)$, which may be difficult to determine due to backgrounds.Comment: 12 pages, LaTeX, no figure

Testing Models of Intrinsic Brightness Variations in Type Ia Supernovae, and their Impact on Measuring Cosmological Parameters

For spectroscopically confirmed Type Ia supernovae we evaluate models of intrinsic brightness variations with detailed data/Monte Carlo comparisons of the dispersion in the following quantities: Hubble-diagram scatter, color difference (B-V-c) between the true B-V color and the fitted color (c) from the SALT-II light curve model, and photometric redshift residual. The data sample includes 251 ugriz light curves from the 3-season Sloan Digital Sky Survey-II, and 191 griz light curves from the Supernova Legacy Survey 3-year data release. We find that the simplest model of a wavelength-independent (coherent) scatter is not adequate, and that to describe the data the intrinsic scatter model must have wavelength-dependent variations. We use Monte Carlo simulations to examine the standard approach of adding a coherent scatter term in quadrature to the distance-modulus uncertainty in order to bring the reduced chi2 to unity when fitting a Hubble diagram. If the light curve fits include model uncertainties with the correct wavelength dependence of the scatter, we find that the bias on the dark energy equation of state parameter $w$ is negligible. However, incorrect model uncertainties can lead to a significant bias on the distance moduli, with up to ~0.05 mag redshift-dependent variation. For the recent SNLS3 cosmology results we estimate that this effect introduces an additional systematic uncertainty on $w$ of ~0.02, well below the total uncertainty. However, this uncertainty depends on the samples used, and thus this small $w$-uncertainty is not guaranteed in future cosmology results.Comment: accepted by Ap

Exploring CP Violation with B_d -> D K_s Decays

We (re)examine CP violation in the decays B_d -> D K_s, where D represents D^0, D(bar), or one of their excited states. The quantity $\sin^2(2\beta + \gamma)$ can be extracted from the time-dependent rates for $B_d(t) -> {\bar D}^{**0} K_s$ and $B_d(t) -> D^{**0} K_s$, where the $D^{**0}$ decays to $D^{(*)+}\pi^-$. If one considers a non-CP-eigenstate hadronic final state to which both D(bar) and D^0 can decay (e.g. $K^+\pi^-$), then one can obtain two of the angles of the unitarity triangle from measurements of the time-dependent rates for $B_d(t) -> (K^+\pi^-)_{D K_s}$ and $B_d(t) -> (K^-\pi^+)_{D K_s}$. There are no penguin contributions to these decays, so all measurements are theoretically clean.Comment: 15 pages, LaTeX, no figure

Type II-P Supernovae from the SDSS-II Supernova Survey and the Standardized Candle Method

We apply the Standardized Candle Method (SCM) for Type II Plateau supernovae (SNe II-P), which relates the velocity of the ejecta of a SN to its luminosity during the plateau, to 15 SNe II-P discovered over the three season run of the Sloan Digital Sky Survey - II Supernova Survey. The redshifts of these SNe - 0.027 < z < 0.144 - cover a range hitherto sparsely sampled in the literature; in particular, our SNe II-P sample contains nearly as many SNe in the Hubble flow (z > 0.01) as all of the current literature on the SCM combined. We find that the SDSS SNe have a very small intrinsic I-band dispersion (0.22 mag), which can be attributed to selection effects. When the SCM is applied to the combined SDSS-plus-literature set of SNe II-P, the dispersion increases to 0.29 mag, larger than the scatter for either set of SNe separately. We show that the standardization cannot be further improved by eliminating SNe with positive plateau decline rates, as proposed in Poznanski et al. (2009). We thoroughly examine all potential systematic effects and conclude that for the SCM to be useful for cosmology, the methods currently used to determine the Fe II velocity at day 50 must be improved, and spectral templates able to encompass the intrinsic variations of Type II-P SNe will be needed.Comment: Accepted for publication by ApJ; data used in this paper can be downloaded from http://sdssdp47.fnal.gov/sdsssn/photometry/SNIIp.tgz; citation errors correcte

Type Ia Supernova Properties as a Function of the Distance to the Host Galaxy in the SDSS-II SN Survey

We use type-Ia supernovae (SNe Ia) discovered by the SDSS-II SN Survey to search for dependencies between SN Ia properties and the projected distance to the host galaxy center, using the distance as a proxy for local galaxy properties (local star-formation rate, local metallicity, etc.). The sample consists of almost 200 spectroscopically or photometrically confirmed SNe Ia at redshifts below 0.25. The sample is split into two groups depending on the morphology of the host galaxy. We fit light-curves using both MLCS2k2 and SALT2, and determine color (AV, c) and light-curve shape (delta, x1) parameters for each SN Ia, as well as its residual in the Hubble diagram. We then correlate these parameters with both the physical and the normalized distances to the center of the host galaxy and look for trends in the mean values and scatters of these parameters with increasing distance. The most significant (at the 4-sigma level) finding is that the average fitted AV from MLCS2k2 and c from SALT2 decrease with the projected distance for SNe Ia in spiral galaxies. We also find indications that SNe in elliptical galaxies tend to have narrower light-curves if they explode at larger distances, although this may be due to selection effects in our sample. We do not find strong correlations between the residuals of the distance moduli with respect to the Hubble flow and the galactocentric distances, which indicates a limited correlation between SN magnitudes after standardization and local host metallicity.Comment: Accepted for publication in The Astrophysical Journal (33 pages, 5 figures, 8 tables